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ABSTRACT 


The fins collector design for solar thermal has widely been used and it has a 
higher thermal efficiency than without Fins. Photovoltaic thermal (PY/T) 
system produced Electrical and thermal energy instantaneously. 
Mathematical modeling based on steady-state thermal analysis of PV/T 
system with and without fins was conducted with matrix inversion method. 
The value results show that the PV/T system with fins collector is higher 
thermal and electrical efficiency than without fins. 
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1. INTRODUCTION 

The collector design is greatest vital units of solar collector system or PV/T system to collect heat 
from solar energy. The use of fins absorber has been conducted widely in the solar collector. Solar collector 
with two-pass fins have been conducted by Naphon[l] with thermal efficiency is about 30-60%. Fudholi et 
al. [2] have compared solar collector two-channel with and without fins performance and price welfares 
analysis. The energy efficiency and cost-effective of two-channel with fins are higher than without 
fins collector. 

Taghi et al. [3] have conducted the analysis of heat reaction for point covered topology in the 
photovoltaic thermal system. The basic idea is the exploit electrical energy through a new point covered 
topology. The experimental investigation shows that the performance of electrical analysis increases about 10 
degrees. Kumari and Babu [4]have compared between theoretical and simulation approach of the 
photovoltaic cell by Matlab-Simulink Situation. The main of purpose is to catch factor of the nonlinear I-V 
calculation by modifying the curve at three ideas: open circuit, extreme power, and short circuit. Sharma et 
al. [5]haveanalyzed the demonstrating and simulation approach of off-grid power generation system by 
photovoltaic. The objective this research is to confirm steady process during responsibility and numerous 
network instabilities in network and islanding connected mode. Zohri et al. [6] have analyzed by 
mathematical modeling of photovoltaic thermal (PVT) with and without a v-groove collector, the energy 
analysis result displays that performance energy with v-groove is higher than without v-groove collector. 

The PV/T system has well in the future because of generating both thermal and electrical energy 
which it’s higher dependability and lower ecological impact [7].Thermal application and different types of 
PV/T system solar collector on basis of physical assortment, scheme condition, and different heat removal 
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improvement technologies have been reported by Tian et al. [8]. Two modes PV/T air heating and PV/T 
water heating of tri-functional PV/T collector has been designed by Guo et al. [9] .The comparative results 
between the simple photovoltaic cell and unit model have been conducted by Zaoui et al. [10] on the 
elementary of connection temperature and isolation. The greater exergy performance with different types of 
PV/T system has been compared to another system[ll]. The special effects of involuntary convection on 
solar cell temperature have been presented by experimental investigation. The solar cell temperature is 
powerfully prejudiced by the competence of drying [12]. 

The electrical and thermal performances of PV/T system based on air collector have been projected 
by Sarhaddi et al. [13]. In evolving the hybrid PV/T fleeting model, the experimental data and theoretical 
modeling have a good arrangement by typical steady-state current-voltage specific curve [14]. The effects of 
dissimilar factors like; length, channel penetration, fluid flow rate and packing factor on electrical and 
thermal performance have been developed by Elsafi et al. [15]. The Combination of thermal collector and 
photovoltaic panel are called PV/T system. A review work has been reported by Good et al. [16]. Using 
glazed PV/T system has been conducted for improving the whole exergy performance by Genetic Algorithm- 
Fuzzy system method [17]. Mohammad et al. [18] have conducted improved the model of solar photovoltaic 
(PV) array along with the implementation of fuzzy logic as maximum power point tracking (MPPT). Surya 
and Sai [19] have done mathematical analysis and simulation of photovoltaic cell using Matlab-Simulink. 
The model analysis based on circuit equation of photovoltaic solar cells with solar radiation and temperature 
parameters. 

Mathematical modeling is very urgent to predict parameters, outlet temperature, thermal and 
electrical efficiency before doing experimental investigation 


2. MATHEMATICAL MODEL 

Based on energy balance, Figure 1 shows schematic heat transfer coefficient for PV/T system with 
fins collector. The structure of collector with and without fins is same principally. For PV panel size (1.2 m x 
0.53 m), fins size 1.2 m of length and 0.03 m of width. The number of fins is 50. 
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Figure 1. Heat Transfer Coefficient for PV/T System with Fins Collector 


For PV/T system without fins is used matrix 3 x 3 to calculate Temperature PV panel T pv , 
temperature fluid 7}and Temperature bottom plate T b . by matrix inverse following: 

[A][T] = [C] 

For PV panel: 

ra(l - tj c ')I = U t (T pv - T a ) + h cl (T pv - T f ) + h r {T pv - T b ) (1) 

For air channel: 

2mC{Tf - Ti)/WL = h cl (T pv - 7>) + h c2 (T b - 7)) (2) 

For bottom plate: 

hr(T pv - T b ) = U b (T b - T a ) + h c2 (T b - 7>) (3) 


Photovoltaic-Thermal (PVT) System with and Without Fins Collector: ... (Muhammad Zohri) 














1758 n 


ISSN: 2088-8694 


(U t + h cl + h r ) —h c i 

h c i Qia + h c 2 ) 

h r h c 2 


—h r 

h C 2 

~(h r + ft c2 + £4) 



r t i 


Ut^a "f TJ cell)I 


i pv 


2 mC 


T r 

T 

— 

WL )Ti 

- 

L *b J 


-U b T a J 


The PV/T system with fins collector is used matrix 3 x 3 to calculate Temperature PV panel T pv , temperature 
fluid Tf and Temperature bottom plate T b . by matrix inverse following: 

[A][T] = [C] 


For PV panel: 

ra(l - r] cell )I = U t (T pv - T a ) + h cl (T pv - 7}) + h r (j pv — T b ) + Q n (4) 

For air channel: 

2mC(T f - T t )/WL = h cl (T pv - 7>) + h c2 (T b - 7>) + Q n (5) 

For bottom plate: 

h r ,(T pv - T b ) = U b (T b - T a ) + h c2 (j b - T f ) (6) 
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Where, 

For Efficiency of fins, collector is 

Qn ~ ~7~ (2kA n lh c ) 0 S tanMH(T f - T t ) (7) 

A fin 

For Electrical efficiency [14] 

T] pv — T]ref \^ — Prefi^pv ~ 4e/)] (^) 

For the temperature coefficient p re f can be written as [20] 

Pref = (T 0 -T ref ) (9) 

For Thermal efficiency by following [21]. 

rith = rhC(T 0 - T t )/IA (10) 


For The convective heat transfer coefficients are given as [22]: 


which, 


h = ——Nu 

D h 


n _ 4 Wd 

h ~2 (w+d) 


( 11 ) 


( 12 ) 


Where, W ,d , D h are the width, high, equivalence diameter of the channel, k is air thermal conductivity and 

Nu. is a Nusselt number. Nusselt numbers are given as, for Re<2300 (laminar flow region): 
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channel, k is air thermal conductivity and Nu is a Nusselt number. Nusselt numbers are given as, for 
Re<2300 (laminar flow region): 
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For 2300 < Re < 6000 (transition flow region): 
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For Re >6000 (turbulent flow region): 

Nu — 0.01 8Re° 8 Pr 04 

Where Re and Pr are the Reynolds and Prandtl number given as: 


For 2300 < Re < 6000 (transition flow region): 


Nu = 0.11 ((Re 3, -12 5)Pr 
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WhereTfe and Pr are the Reynolds and Prandtl number given as: 


For Re > 6000 (turbulent flow region): 
Nu = 0.01 8Rc" 8 Pr 04 


Re= mA 


Pr= 


Aa 

juC 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 


The theoretical model simulated that for a short collector or less of 1 m. Then, the mean air temperature is 
then equal to the arithmetic mean, where: 

T fo+rj 

' 2 ( 20 ) 

The physical properties of air are hypothetical vary linearly with temperature ( °C) by Fudholi [23]: 

Specific heat 

C = 1.0057 + 0.000066(7 — 27) (21) 

Density, 

p = 1.1774 - 0.00359 (T - 27) (22) 

Thermal conductivity 

k = 0.02624 + 0.0000758 (T - 27) (23) 

Viscosity 

p = [1.983 + 0.00184(7 - 27)]10" 5 (24) 

The heat transfer coefficient to wind according to Ong [24] is 

h w = 2.8 + 3.3V (25) 

Where, h w heat transfer coefficient due to wind and V is the wind velocity. The heat transfer coefficient from 
panel cell to sky 

t _ crgpy ( Tp V +T s ) ( Tpy+Ts ) (Tp V ~T s ) 


, t T> p VS 

l pv~ l a 
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Where T s is the sky temperature, T c is the photovoltaic panel temperature. 
T s = 0.05527 a 15 


(26) 

(27) 

(28) 


3. RESULTS AND ANALYSIS 

Figure 2 shows mass flow rate versus outlet temperature of air at PV/T system with and without fins 
collector with solar intensity of 600 W/m 2 and 800 W/m 2 . The result indicates that by growing the mass flow 
rate simultaneously fallen the temperature outlet of air. The mass flow rate from 0.01 kg/s to 0.02 kg/s is 
laminar flow state with (Re <2300) and mass flow rate from 0.03 kg/s to 0.05 kg/s is transition flow state 
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with (2300 < Re < 6000). The maximum temperature outlet at intensity of 800 W/m 2 with fins collector is 
41.39 °C and the minimum temperature outlet at intensity of 600 W/m 2 is 29.77°C. 


43.00 



-•-With Fins (600 W/m2) Without Fins (600 W/m2) 

—A— With Fins (800 W/m2) —Without Fins (800 W/m2) 


Figure 2. The Mass Flow Rate Versus Outlet 
Temperature of With and Without Fins Collector 
at Solar Intensity of 600 W/m 2 and 800 W/m 2 



—♦—With Fin (600 W/m2) -B-Without Fins (600 W/m2) 
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Figure 3. The Mass Flow Rate Versus Thermal 
Efficiency of With and Without Fins Collector 
at Solar Intensity of 800 W/m 2 and 600 W/m 2 



Figure 4. The Mass Flow Rate versus Cell Efficiency of Fins Collector 
for the Solar Intensity of 600 W/m 2 and 800 W/m 2 


Figure 3 shows the mass flow rate versus thermal efficiency of with and without fins collector. For 
the mass flow rate from 0.01 to 0.05 kg/s, the thermal efficiency maximum at the solar intensity of 600 W/m 2 
and 800 W/m 2 is 43 % with fins collector and thermal efficiency minimum at the intensity of 600 W/m 2 and 
the intensity of 800 W/m 2 is 26% without fins collector. When mass flow rate value decreases thermal 
efficiency follow it. This situation can be concluded that the higher the mass flow rate, the higher thermal 
efficiency. Figure 4 shows distribution mass flow rate versus electrical efficiency with and without fins 
collector. The conduct of electrical efficiency indicates that electrical performance increase with following 
increasing mass flow rate. The use of fins collector be able to upsurge the electrical efficiency. It show that 
the effect of fins collector make cooling photovoltaic panel. 


Table 1. The Comparison with the Other Collector Designs in Reference [25,26] 


Designs of collector 

Predicted Results 

Thermal efficiency Electrical efficiency 

glass 

36% 

- 

metal sheet collector and glass 

41% 

- 

fins and glass 

51% 

- 

Two-absorber (no cover type ) 

65% 

8.4% 

Two-absorber (cover type) 

66% 

8.5% 

Free channel 

64% 

8.6% 

current study (fins collector) 

43% 

14% 
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Table 1 demonstrates the comparison designs of collector between using fins collector, metal sheet, 
with glass, two-absorber with and without cover type. The usage of fins is higher thermal efficiency than 
without fins and with metal sheet. The other way, PVT with fins collector combination with glass cover and 
two-absorber with and without cover type are higher thermal efficiency than just with fins collector without 
glass cover. It shows that the glass cover and two-absorber be able to rise drying up in PVT system. 

4. CONCLUSION 

Mathematical modeling for predicting thermal and electrical efficiency of photovoltaic thermal 
(PV/T) system with and without fins has been conducted. The maximum results of outlet temperature, 
thermal efficiency, and electrical efficiency are 41.39 °C, 43%, 14% respectively. Using fins collector to 
cooling PV panel has higher efficiency than without fins collector. 


5. 

APPENDIX 




Nomenclature 

A 

area 

m 2 

C 

specific heat of air 

J/kg.°C 

d 

channel high 

m 

h 

heat transfer coefficient 

W/m 2 .°C 

L 

length collector 

m 

I 

intensity 

W/m 2 

w 

width collector 

m 

Pr 

Prandtl number 


Re 

Reynold number 


T 

Temperature 

°C 

Greek letters 


£ 

emissivity 


t 

transmission coefficient 


a 

absorption coefficient 


P 

dynamic viscosity 


V 

efficiency 


Subscripts 


i 

inlet 


0 

outlet 


f 

fluid 


s 

sky 


r 

radiation 


c 

convection 


b 

back plate 


a 

ambient 


pv 

photovoltaic panel 


a 

ambient 
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